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Abstract 

Electric charge and color breaking minima along third generation squark 
directions do appear in the Minimal Supersymmetric Standard Model for 
particular regions of the corresponding supersymmetric parameters At, fi, 
tan/3, rriQ and mfj. We have studied possible instabilities of electroweak 
sphalerons along non-trivial squark configurations. We have found that in- 
stabilities along the latter lie in the region of parameter space where the 
electroweak minimum is not the global one. Thus charge and color con- 
servation imply that the standard electroweak sphaleron is not destabilized 
along squark directions. 
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1. The existence of sphalerons (static and unstable solutions to classical field equations) 
in the SU(2) gauge theory with Higgs fields in the fundamental representation has been 
known since long ago , triggering the hope of generating the baryon asymmetry of the 
Universe at the electroweak phase transition 0. The value of the sphaleron energy, 
an essential ingredient for baryogenesis mechanisms 0, ^, ^ has been the object of 
detailed numerical calculations in the Standard Model ||, ^, |10|, and translates into an 
upper bound on the Higgs boson mass well below the experimental lower limit. This 
bound comes from the weakness of the phase transition in the Standard Model. 

The minimal supersymmetric extension of the Standard Model (MSSM) [|lT| is a 
very appealing candidate, both from the theoretical Q and the experimental side, to 
describe physics at the electroweak scale. The strength of the phase transition in the 



MSSM has been extensively studied in the hterature |[T^, |T^, |T^. Recently it has been 
proven that the phase transition can be considerably strengthened in the MSSM in 
the presence of light supersymmetric partners of the right-handed top quark (stops) 



and small values of tan/5 [^, 16, 17|, opening thus the window for baryogenesis at the 



electroweak phase transition. These results have been recently confirmed by an explicit 
calculation of the sphaleron solution and the sphaleron energy in the MSSM, at zero 
and finite temperature, including all relevant one-loop radiative corrections . 



The sphaleron solution presented in Ref. is the natural generalization of the 
Standard Model sphaleron solution and involves the two Higgs doublets and the SU(2) 
gauge fields. However, the main feature of supersymmetric models is the existence of a 
plethora of new scalar fields, the sfermions, which might affect the sphaleron solution 
if non-trivial configurations turn out to be energetically favoured. In particular, for the 
region of supersymmetric parameters where the phase transition is strengthened, there 
are electric charge and color breaking (CCB) local minima along the left-handed 
and/or right-handed stop directions, which might destabilize the standard sphaleron in 
the MSSM and then endanger previous results based on standard sphalerons and the 
analysis of the phase transition. 

In this letter will study the effect of possible non-trivial sfermion configurations on 
the sphaleron solution and its energy. We find that as far as electric charge and color 
are not broken along squark directions, i.e. as far as the standard electroweak minimum 
along the neutral component of the Higgs fields remains as the global minimum, the 
trivial (zero) configuration for sfermion fields is always energetically favoured. 



2. Neglecting sfermion mixing between different generations (as suggested by con- 
straints from flavour changing neutral current processes), we can simplify the problem 
by considering just one generation of squarks and sleptons. On the other hand, as a first 
approximation, we will fix to zero the slepton fields. The consistency of this fixing is 
protected by a global symmetry and will be justified a posteriori by our results. In fact 
we will work with third generation squarks, whose large Yukawa couplings would allow 
tunneling from the symmetric vacuum at temperatures above the electroweak phase 
transition temperature, or from the standard non-symmetric electroweak vacuum, for 

^Yielding a 'technical' solution to the hierarchy problem. 

^LEP precision measurements leading to gauge coupling unification. 
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temperatures below, to a possible lower CCB vacuum. We will also work in the approx- 
imation of taking 5(1 = so that the U{1)y gauge field can be consistently set to 
zero. This allows a spherically symmetric ansatz. So, we are looking for a sphaleron-like 
solution described by 



and 
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The lagrangian density for these fields is given by 

^ = -lW^,W-'''^ + {D,H,y{D^^H,) + {D^H2y{D^^H2) + {D^QL)HD''QL) 
+ {d,UI,y {d>^U%) + {d,Dl)^ (d^D'k) - V,^{H,, H2, Ql, U'k, D%) 



(3) 



where VF^j, stands for the SU{2) field strength and is the effective potential. We 
can expand the effective potential (at zero temperature) as 



Kff = Vo{H,, H2, Ql, u'r, oy) + V,{H,, H2, Ql, U'r, 0%) + 



(4) 



where Vq is the tree-level potential and Vi contains the one-loop radiative corrections. 
The tree-level potential can be decomposed as: 



with 
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the SU(2) and SU(3) D-terms. Vp is determined from the superpotential 
W = htQL ■ H2U'r + hQL ■ H.D'j, + fiHi ■ H2, 



and reads: 



IhQLD'j.-fiH^l^ + \htQLU^^ + fiH, 



The soft-breaking terms are: 

Voft = m\ hIHi + ml HIH2 + ml (Hi ■ H2 + h.c.) 

+ nil QIQl + U'JUl, + ml DiD'i, 

+ {htAt Ql ■ H2Uf, + hA, Ql ■ H^D% + h.c.) 



(9) 



(10) 



where m\, m\ and m|, can be traded in favour of the supersymmetric parameters, 
tan/3 and m^-, using the minimization conditions of the radiatively corrected effective 
potential. Finally, for the one-loop effective potential one should use its expression in 
the DR renormalization scheme: 



^1 
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(11) 



where Q is the renormalization scale. We have taken the approximation where only 
Higgs background fields are considered in the one- loop correction, Vi = Vi{Hi,H2). 
This correction can be absorbed in a redefinition of the tree-level parameters in the 
Higgs sector |^ and hence provides large corrections to Higgs boson masses. On the 
other hand the tree-level potential along the squark directions is dominated by the 
strong and top-quark ht Yukawa couplings. Since the minima of the potential will 
always be located at field values of the order of the weak scale, we do not expect large 
radiative corrections for them and we can safely put to zero the background fields Ql, 
U^^ and in (|ll|). 

Choosing the temporal, radial gauge we can write the static, spherically symmetric 
ansatz 
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for the sphaleron 0, where = + + and all functions are real. Notice that we 
have not included gluons, while Ql and U'^ are aligned on a unique SU(3) direction, 
which provides a self consistent ansatz. 

Replacing the ansatz (0) in @ we obtain in a straightforward fashion the scalar 
potential as a function of the background fields hi, h2, q and u: 

V^s = Ves(hi,h2,q,u). (13) 

We have explored the region where the supersymmetric parameters 

{rriQ, mu, rriA, tan (3, At, fi) 

can take all values allowed by experimental and theoretical constraints. In particular 
large values of mq are preferred by LEP precision measurements, while large values of 
niA, triA ~ TTiQ, small values of tan/5, tan/5 < 3, and moderately negative values of m^, 
mlj > —mf, are favoured by the strength of the electroweak phase transition, and will be 
focused on in the following discussion. This region has been recently proven to appear 
naturally from the usual radiative scenario of electroweak symmetry breaking ||2^ . On 



the other hand it is clear that the presence of CCB minima and hence of possible 
instabilities of the standard sphaleron along the squark configurations will be highly 
favoured by mfj < values. 

The analysis of the potential ([l3|) reveals the presence of the standard electroweak 
minimum, at q = u = 0, hi = v cos /?, /i2 = f sin j3, where v = 174.1 GeV is the standard 
Higgs vacuum expectation value (VEV), with a depth 

= -\mlv'' (14) 



where niH is the lightest Higgs boson mass in the limit niA oo [|18|. For the case we 
are considering of large values of rriA, ran is just the lightest Higgs boson mass. 
For values of At = At + /x/ tan/? < At max, with 

max _ _ gS^riH ^^^^ 

\/3ht sin PiriH 



m, 



2 



there is a minimum at hi = h2 = q = 0, = — where mfj = —mfj. The condition 

93 _ 
that this CCB minimum not be lower than the standard minimum imposes on niu the 



upper bound niu < m?[^* given by fTS 



s«=(<|My'* (16) 

We have plotted in Fig. 1 (long-dashed line) m^'* given by Eq. (0) as a function of At 
for rriQ = rriA = 500 GeV and tan (3 = 3. We can see the curve ends at a value of At 
given by Eq. (|T3p. 

^Notice that we have fixed the right-handed squark fields to the trivial configuration Dl — -D^ = 0. 
This result is consistent with the smallness of the Yukawa coupling hb, that we will neglect from here 
on, and will be justified a posteriori by the numerical results of this work. 
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For At > At max there is a CCB minimum with all fields hi, h2, q and u acquiring 
VEVs. The location of the minimum can only be computed numerically and m^^*, 
at which the CCB minimum becomes degenerate with the standard electroweak mini- 
mum is shown in Fig. 1 (short-dashed line) for the same values of the supersymmetric 
parameters. 

In short, in the region below (above) the dashed line the electroweak minimum (the 
CCB minimum) is the global one, and only below the dashed line is electric charge and 
color guaranteed to be unbroken. For illustration we have depicted also in Fig. 1 the 
(thin solid) line where the right-handed squarks are massless. Then only below the thin 
solid line there exists a (local or global) standard electroweak minimum. 



The last issue we want to address is the possible (in) stability of the standard 

It is again clear that the region < 0, where the 
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sphalerons in the MSSM |18 



right-handed squark direction is unstable at the origin, favours the presence of such 
instabilities. One should now substitute the ansatz (|12]) into the Euler- Lagrange equa- 
tions, 

■92 



■ yz \ - 

^ Y=Hi,2,Ql 



dV 



Y^(T^{D^Y) - (L)'^F)V'^F 



(17) 



9Xt 



where X = H12, Ql, Uf^, supplied by the appropriate boundary conditions. In some 
cases we expect to end up with a solution described by vanishing squark fields, whereas 
in other situations non-vanishing profiles can be preferred. Before embarking ourselves 
in the complete solution of the equations of motion, we will study the possible standard 
sphaleron instabilities along the squark field configurations. Hence we will treat the 
standard (squarkless) sphaleron as a background and we will study the corresponding 
squark equations in this background. This approximation will be justified a posteriori, 
and will provide us with an intuitive picture of the parameter space region where squarks 
are involved in the sphaleron solution. 

The energy of the squark fields in the sphaleron background is given by: 



Att / dr {drqf + ^(1 - ffq^ + {drU 



+ 



^q' + 2q\hl-hl)\+^-l-{q' - u'f 

-,272—2 , L2r2~2 



+ h^q ^ + 2fihtq u hi + hth2U ^ + h^h^q 
+ "^Q^ ^ + 'm^iju'^ + 2htAtq h2u' 

Non-vanishing squark field configurations will be associated to negative values of 
the previous functional. We can split it as follows: 



(19) 
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where the first piece includes quadratic terms in the squark fields and the second one 
collects the remaining, quartic, ones. Notice that S^E~ > and then a decreasing of 
the energy is possible only in the case where the quadratic terms are negative. These 
are given by: 



where 



5'E~=47r Idr 



rq{r) ru{r) 



rq{r) 
ru{r) 



(20) 



-9r + - ff + |(/^? - hi) + hfhl + ml h^i^Th + A,h„ 



ht{fihi + Ath2 



-dl + hlhl + ml 



(21) 



To determine the instabilities of the standard sphaleron along the squark direction 
one must diagonahze the quadratic operator in (^). Each eigenvector of with a 
negative eigenvalue 

\ 1 = -^4 1 (22) 

ru{r) ru[r) ^ ' 

generates an instability. 

We have studied the spectrum of the operator ( PT| ) numerically, using a discretized 
representation of the quadratic energy functional. The MSSM sphaleron is characterized 
by just one scale ~ My/-, since it is mainly controlled by the W mass and the lightest 
Higgs boson mass, the latter never being hierarchically larger than the former Q. When 
squarks are included, a second scale, M, could emerge in our solution. It would be 
roughly given by the inverse radius of the region where squarks are non- vanishing. The 
value of this scale depends on the spectrum of our MSSM model, but we expect it to be 
in the range Mw ^ M < Msusy, where Msusy is typically ^ 1 TeV. In order to cover 
both scales, we will follow the approach of Ref. [1^ and use a uniform discretization of 
the variable: 

1 + Mr 



.1 + Mwr] 

where s goes from to 1 when r varies from to infinity. We have done our analysis 
using 200 points and taking M values within the above mentioned range. We have 
verified, by increasing the number of points and comparing, that the lower eigenvalue 
is calculated with and error smaller than 0.001 in units. 

We have plotted in Fig. 1, for the same values of the supersymmetric parameters 
(tan/3, m^, m^) as before the (thick solid) line along which there appears a negative 
eigenvalue along a non-trivial squark configuration. So above (below) the thick solid 
line the standard sphaleron in the MSSM is unstable (stable). We can see from Fig. 1 
that the instability appears only in the region where the CCB minimum is the global 
one and electric charge and color are broken. We have stopped the instability line when 
it intersects the thin solid line, where the standard electroweak minimum disappears. 



^This situation is different in the Standard Model, where the quartic couphng of the Higgs field is, 
in principle, arbitrary, so that it can be arbitrarily heavier than the W boson, in which case a second 
scale, Mh, appears. 
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We have not shown the region with mfj > 0, but it becomes clear from the behaviour 
of the thick sohd hne in Fig. 1 that there is no instabihty region for mfj > and any 
value of Af. 

We have also explored general values of the supersymmetric parameters in the region 
1 < tan P ^ 20 (where we can comfortably neglect the bottom Yukawa coupling hb as 
compared to ht), 100 GeV < tua ~ ttiq and 200 GeV < rriQ < 500 GeV, and we have 
found a qualitative behaviour similar to that in Fig. 1. Even in the region of large 
values of tan /3, where hb can no longer be neglected, our results holds since the possible 
instability along Dl, D"^ would be controlled by the entry 

-d'^ + hlhl + ml 

in the corresponding quadratic operator, similar to that in (pTD , which in turns depends 
on the value of the bottom-quark mass. Therefore, the present experimental bounds 
on supersymmetric masses impose mj^ > which prevents, using the results obtained 
in this paper, the existence of instabilities along the squark fields Dl,D1^. Since the 
previous results obviously apply to slepton fields, our choice of the ansatz where slepton 
and down-squark fields vanish remains fully justified. 



4. Using the previous results, which concern the MSSM standard electroweak sphaleron 
at zero temperature, we can infer the stability of the standard sphaleron against non- 
trivial squark configurations at temperatures higher than, or equal to, the critical tem- 
perature of the electroweak phase transition, Tc ~ 100 GeV. The argument goes as 
follows: The main effect of finite temperature corrections in the MSSM sphaleron en- 
ergy and profiles can be encoded, to a reasonable accuracy, in the (asymptotic) vacuum 
expectation values of the Higgs fields, i.e. vi(T) and V2iT) |jl8[. On the other hand, the 



only modification in the quadratic operator ( ^iD from finite temperature corrections is 
to replace the squark masses mg, mfj by effective thermal masses, mqiT) = mQ+IlLiT) 
and m^(T) = mfj + Ilji{T), where the self-energies IIl{T) and IIfi{T) can be found in 
Ref. Ijl^, and whose precise expressions are not important for the sake of the present 
discussion. A necessary condition at T > Tc is that mqiT) > 0, m\j{T) > since 
otherwise the phase transition would proceed at some temperature Tccb > Tc towards 



the finite temperature CCB minimum through a strong first order phase transition [15 



Using the results of the present work we can infer there would be no instability region 
for the range of temperatures T > T^, since negative squared masses have been found to 
be a necessary condition for the existence of instabilities. In other words, the stability 
of the standard electroweak minimum at the critical temperature implies stability of 
the standard electroweak sphaleron along squarks configurations. Of course, below Tc, 
m|(T) and/or m^(T) can become negative and an instability region will start to grow 
up and evolve towards the zero temperature result shown in Fig. 1. 



5. In conclusion, we have studied possible instabilities of the standard electroweak 
sphaleron in the MSSM along the most dangerous third generation sfermion configura- 
tions, with CCB minima, which might decrease the sphaleron energy. At zero temper- 
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ature, we have found that instabihties he in the region of supersymmetric parameters 
where CCB minima are the global minima. Thus imposing electric charge and color 
conservation is a sufficient condition to ensure the stability of the standard sphaleron 
along third generation sfermion directions. Finally, at the electroweak phase transition 
temperature, imposing the stability of the electroweak minimum prevents the appear- 
ance of instabilities along sfermions configurations. 
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Figure 1: For tan/3 = 3 and uiq — tua — 500 GeV: a) Above the thin sohd hne the 

elcctroweak minimum disappears; b) Above (below) the dashed hnes electric charge 
and color are broken (conserved); c) Above (below) the thick solid line the standard 
sphaleron is unstable (stable) along the squark directions. 
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